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cn ■ Abstract 
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In a recent work we described the form taken by the Minkowski-space gluon propagator in the 
presence of an < A^Aa > condensate. This condensate has been seen to play an important role 
in calculations which make use of the operator product expansion. The latter work has led to 
the publication of a large number of papers which discuss how the < A^Aa > condensate could 
play a role in a gauge-invariant formulation. In the present work we discuss the properties of both 
the Euclidean-space and Minkowski-space gluon propagator. In the case of the Euclidean-space 
propagator we can make contact with the results of QCD lattice calculations of the propagator in 
the Landau gauge. With an appropriate choice of normalization constants, we present a unified 
representation of the gluon propagator that describes both the Minkowski-space and Euclidean- 
space dynamics in which the < A^Aa > condensate plays an important role. 
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I. INTRODUCTION 



Recently, studies making use of the operator product expansion (OPE) have provided 
evidence for the importance of the condensate < A a ^A^ > [1-3]. (There is a suggestion that 
such a condensate may be related to the presence of instantons in the vacuum [4].) The 
importance of that condensate raises the question of gauge invariance and there are now a 
large number of papers that address that and related issues [5-19]. We will not attempt to 
review that large body of literature, but will consider how the presence of an < A a ^A^ > 
condensate affects the form of the gluon propagator. We may mention the work of Kondo 
[7] who was responsible for introducing a BRST-invariant condensate of dimension two, 

Q = fl < J d4xTT (\M x )M x ) - aic ( x ) ■ c(s)) >, (1-1) 

where c(x) and c(x) are Faddeev-Popov ghosts, a is the gauge-fixing parameter and Q is 
the integration volume. Kondo points out that Q reduces to A 2 min in the Landau gauge, 
a — 0. The minimum value of the integrated squared potential is A 2 miw which has a definite 
physical meaning [7]. 

In a recent work we considered the Minkowski-space gluon propagator in the presence of 
an A 2 condensate [20]. We found that the propagator has no on- mass-shell poles, so that 
the gluon was a nonpropagating mode in the presence of the vacuum condensate [21]. The 
form we obtained for the propagator was 

(k^k v \ 
r - D(k) (1.2) 

with 
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Here Z 1 is a normalization parameter which we put equal to 3.82 so that we may obtain a 
continuous representation as we pass from Minkowski to Euclidean space. In Fig. 1 we show 
D(k) with m 2 = 0.25 GeV 2 . (We remark that D(k) = when k 2 = m 2 , D(k) = -Z 1 /m 2 
at k 2 = 0, and D(k) — > Z\jk 2 for large k 2 .) If we choce Z\ = 15.28m 2 = 3.82 our result 
for the propagator will be continuous at k 2 = when we consider both the Euclidean-space 
and Minkowski-space propagators. 

The organization of our work is as follows. In Section II we discuss the Euclidean-space 
gluon propagator as obtained in a lattice simulation of QCD in the Landau gauge and in 
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the absence of quark degrees of freedom [22]. (We also record in the Appendix a number 
of semi-phenomenological analytic forms which are meant to represent the Euclidean-space 
propagator.) In Section III we summarize the results of our analysis and provide some 
additional discussion. 



II. QCD LATTICE CALCULATIONS AND PHENOMENOLOGICAL FORMS 
FOR THE EUCLIDEAN-SPACE GLUON PROPAGATOR 



Results for the gluon propagator obtained in a lattice simulation of QCD are given in Ref. 
[22]. In that work the authors also record several phenomenological forms. We reproduce 
these forms in the Appendix for ease of reference. Of these various forms we will make use 
of model A of Ref. [22] which has the form 

AM 2a 1 



with 



D L (k 2 ) 



L(k 2 ,M) 



+ 



(k 2 + M 2 y+ a k 2 + M 2 



-\n(k 2 + M 2 )(k~ 2 + M~ 2 ) 
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-do 



(2.1) 



(2.2) 



and cLd = 13/22. The parameters used in Ref. [22] to provide a very good fit to the QCD 
lattice data are 



2.01 



+4 
-5; 



and 



A = 9.841J8, 



M = 0.54^, 



a = 2.17^9- 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



Note that M in GeV units is 1.018 GeV. Rather than work with the lattice data we will 
use Eqs. (2.1)-(2.6) when we compare our results with the lattice data. In Fig. 2 we show 
k 2 D L (k) of Eq. (2.1) and in Fig. 3 we show D L (k). These functions are represented by the 
solid lines in Figs. 2 and 3. Note that Eq. (1.3) may be written in Euclidean space as 



This form is useful for k\ < 1 GeV 2 and we therefore consider various phenomenological 
forms which may be used to extend Eq. (2.7) so that we may attempt to fit the lattice result 
over a broader momentum range. To that end, we make use of Ref. [23]. The authors of 
that work define the Landau gauge gluon propagator as 

< Al(k)K(k') >= V5(k + k')5 ab (V - ^) (2.8) 

with 

(2 \ 2k 

A2 * ) Hk 2 ))~\ (2-9) 

and 7 = —13/22. (We do not ascribe any particular significance to Eq. (2.9). We use Eq. 
(2.9) as a phenomenological form which could be replaced by a form which provides a better 
fit to the data within the context of our model at some future time. We believe Eq. (2.9) is 
useful, since it is a simple matter to remove the first term of that equation and introduce a 
propagator that has the small k 2 behavior of our model.) 

The authors of Ref. [23] introduce two choices for a{k 2 ) of Eq. (2.9). We use their form 
for a 2 (k 2 ): 

Mk 2 ) = — v ^ - xa21 . (2.10) 
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l QCD , 

In their analysis they put re = 0.5314, A QCD = 354 MeV, a(0) = 2.74, a x = 0.0065 and 
a 2 = 2.40. (Here, we have not recorded the uncertainties in these values which are given in 
Table 2 of Ref. [23].) As we proceed, we will change these values somewhat. As a first step 
we remove the first factor in Eq. (2.9) and write 

Z{k 2 ) = Z 2 {a 2 {k 2 ))-\ (2.11) 

We now use ai = 0.0080 and a 2 = 2.10 rather than the values given above. In Fig. 4 we 
show (a 2 (^)) 13 ^ 22 as a function of k, using our modified values of a\ and a 2 . 
We now define 

DEikE) = - 2 2 ))-; m2 . (2 , 2) 

K E + m 3 kl+mi 

The function —k 2 D E (k E ) is shown in Fig. 2 as a dotted line. In this calculation we have 
put Z 2 = 2.11. We find a good representation of the lattice result for k E < 2 GeV, 




1 1 . 1 . 1 . 1 1 1 . 1 

1 2 3 4 5 6 

k 2 (GeV 2 ) 



FIG. 1: The function D(k 2 ) of Eq. (1.3) is shown in Minkowski space. The value for large k 2 is 
given by Z 1 /k 2 with Z 1 = 3.87. Here m = 0.50 GeV. 

In Fig. 3 we compare D E {k,E) with the result of the lattice calculation which is represented 
by the solid line. In Fig. 5 we combine our results in Minkowski and Euclidean space and 
show the values of k 2 D(k 2 ) for both positive and negative k 2 values. For positive k 2 we use 
D(k) of Eq. (1.3) and for negative values of k 2 we use DE{k 2 E ) of Eq. (2.12). Equality of 
these functions at k 2 = implies Z\ = Z2(a(0)) 13 ^ 22 , or Z\ = I.8IZ2. (In our work we have 
used Zi = 3.82 and Z 2 = 2.11. See Eqs. (1.3) and (2.12).) In Fig. 6 we show D(k 2 ) rather 
than k 2 D(k 2 ), which was shown in Fig. 5. 

III. DISCUSSION AND CONCLUSIONS 

In this work we have provided a representation of the gluon propagator in both Euclidean 
and Minkowski space. The Minkowski-space propagator has only complex poles and that 
implies that the gluon is a nonpropagating mode in the QCD vacuum. Our analysis takes 
into account the important condensate < A^A® > which is responsible for mass generation 
for the gluon. Our work has some relation to that of Cornwall [24] who obtained a gluon 
mass of 500 ± 200 MeV in his analysis. Cornwall also suggested that "quark confinement 
arises from a vertex condensate supported by a mass gap." 

In recent work Gracey obtained a pole mass of the gluon of 2.13Aj^ in a two- loop 
renormalization scheme [25]. If we put Ajfg — 250 MeV, the mass obtained at two-loop 
order in Ref. [25] is 532 MeV, which is close to the value of 500 MeV used in the present 
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FIG. 3: The function —DE(k) is shown. The solid line represents the QCD lattice data, while the 
dotted line represents -D E {k) of Eq. (2.12). [See Fig. 2.] 

work. (We remark that in Ref. [21] we obtained a gluon mass of 530 MeV, if we made use 
of Eq. (3.18) of that reference, which includes the effect of including various exchange terms 
in our analysis of the relevant matrix elements.) 
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FIG. 4: The function (a 2 (k)) 13 / 22 is shown. [See Eq. (2.10).] Note that (a 2 (0)) 13 / 22 = 1.81. 
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FIG. 5: For k 2 > the solid line represents k 2 D(k 2 ) with D(k 2 ) given by Eq. (1.3). Here, 
Z 1 = 3.82. For k 2 < we show k 2 D E (k 2 ), where D E {k 2 E ) is given by Eq. (2.12) with Z 2 = 2.11. 

APPENDIX A 



For ease of reference we record various semi-phenomenological forms which are meant to 
represent the Euclidean-space gluon propagator. 
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FIG. 6: Same as Fig. 5 except that D(k 2 ) is shown. 



Gribov [26]: 



^ 2 > = FT^< M) - 



Stingl [27]: 



D (k ) 



Zk? 



k 4 + 2A 2 k 2 + M 4 



L(k 2 ,M). 



Marenzoni et al. [28]: 



D L (k 2 ) 



(k 2 ) 1+a + M 2 ' 



Cornwall I [24]: 



D L (k 2 ) = Z 



[k 2 + M 2 (k 2 )]\n 



k 2 + AM 2 {k 2 ) 
A 2 



where 



M(k 2 ) = M 



In 



fc 2 +4M 2 
A 2 



In(^) 



-6/11 



Cornwall II [29]: 



D L {k 2 ) = Z 



[k + M ] In 



A; 2 + 4M 2 ' 
A 2 



Cornwall III [29]: 



D L (k 2 ' 



k 2 + Ak 2 \n(^) + M 2 ' 



Model A [22]: 



D (k ) = Z 



AM 



2a 



1 



{k 2 + M 2 ) 1+a k 2 + M 2 

The parameters for model A are given in Eqs. (2.3)-(2.6). 
Model B [22]: 

AM 2a 1 



L(k 2 ,M) 



D (k ) = Z 



+ 



_(k 2 ) 1+a + (M 2 ) 1+a k 2 + M 2 



; L(k 2 ,M) 



Model C [22]: 



D L (k 2 ) = Z 



A fl,2/]\ /r 2\a 



M 2 



e -(k*/M*y + 



k 2 + M 2 



L(k 2 ,M) 



(A7) 



(A8) 



(A9) 



(A10) 
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